Asymptotically Flat Initial Data for 
Gravitational Wave Spacetimes, 
Conformal Compactification 
and Conformal Symmetry 



66 ! Sascha Husa 



Department of Physics and Astronomy, University of Pittsburgh 



> ; October 11, 1998 

O 

(N 



Abstract 

We study the utilization of conformal compactification within the con- 
l/-) ■ formal approach to solving the constraints of general relativity for asymp- 

. totically flat initial data. After a general discussion of the framework, 

particular attention is paid to simplifications that arise when restricting to 
a class of initial data which have a certain U(l)xU(l) conformal symmetry. 
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$L 1 Introduction 

i 

kjrjj Applying the conformal method Jj] of solving the constraints of general relativ- 

ity to asymptotically flat (AF) initial data, the usual choice is to pick an AF 
representative of the conformal equivalence class as base metric. Here instead 
we discuss the alternative of choosing a conformally compactified seed geometry. 
This formalism offers a number of advantages, e.g. it is not necessary to introduce 
any artificial cutoff at some finite radius in numerical work, and all asymptoti- 
cally flat regions in multiple-black-hole initial data can be treated on an equal 
footing. Furthermore the method is well suited to exploit a certain conformal 
symmetry, which yields further simplifications. 

The conformal symmetry considered here consists of a U(l) x U(l) group 
of conformal isometries associated with two commuting, orthogonal conformal 
Killing vector fields, and will be called (conformal) toroidal symmetry. Their 
action extends to a U(l) x U(l) symmetry on the many-point compactification 
S 3 of the physical initial hypersurface. 

Based on the discussion presented here, a numerical analysis, including the 
existence and properties of apparent horizons, has been carried out for (i) Brill 
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waves (plus black holes), (ii) initial data containing a marginally outer trapped 
torus, and (iii) time-asymmetric initial data, where the extrinsic curvature is ob- 
tained as an exact solution for non-conformally flat geometries with conformal 
symmetry. These results have been discussed in [H and B. The organization of 
this paper is as follows: first we briefly review the well known formulation of the 
constraints on a compactified background in Sec. [| then we discuss toroidal sym- 
metry and the usual Brill ansatz for axisymmetric initial data in its compactified, 
toroidally symmetric incarnation in Sec. 0. Section ^ discusses the constraints in 
the context of toroidal symmetry, in particular the case of time symmetry, while 
Sec. [| deals with the case of multiple black holes. 



2 Constraints and Compactification 
2.1 Time Symmetry 

In this section, we will briefly discuss the formulation of the constraints on a 
compactified background. The standard 3+1 split of Einstein's equations yields 
one vector (momentum) and one scalar (Hamiltonian) constraint. Choosing a 
maximal slice (defined as a slice of vanishing mean curvature) as the initial Cauchy 
surface E, the Hamiltonian constraint reads 

3 R{hab\ = K ab K ab , (1) 

where 3 R[h a b] is the scalar curvature of the spatial 3-metric h a b, and K a t, is its 
extrinsic curvature. The momentum constraint on a maximal slice is 

D b K b a = 0. (2) 

For simplification the time-symmetric case (K a b = 0) is discussed first. Then 
the momentum constraint (0) is satisfied identically, and only the Hamiltonian 
constraint ([]]) remains. In the conformal approach to solving the constraints, 
the physical metric is constructed from a freely specifiable seed metric h a b via 
conformal rescaling: h a b = ip 4 hab- Assuming h a b to be asymptotically flat this 
ansatz transforms the Hamiltonian constraint into the linear elliptic equation 

Lhi> = (-A s + -i^J i/> = 0, V > 0, V -> 1 at oo, (3) 

for the conformal factor if). The differential operator is called the conformal 
Laplacian of the metric h. 

Note that, since the mass can be given an arbitrary value by conformal rescal- 
ings, we may in particular choose a base metric with zero mass for simplicity. The 
mass is then fully encoded in the conformal factor that solves the Hamiltonian 
constraint. For a massless metric h, there exist regular one-point compactifica- 
tions (£, g) of (£, h). So g a b = oj 2 h a b is a smooth metric on E = S - A, where A 
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is the point "at infinity", it can be extended smoothly to (E,g) and u is a regular 
asymptotic distance function (RADF) of (£,<?) at A (compare e.g. Q and j5|). 
A function uj is called an RADF near a point A of a manifold (S,g) if it is C°°, 
positive and satisfies 

cu| A = 0, D a uj\ A = 0, {D a D b u - 2g ab )\ A = 0, D a D b D c u\ A = 0. (4) 

In asymptotically Cartesian coordinates on E, the asymptotic behavior ofu; is 
given by u = 0(r~ 2 ) near infinity. 

Going over to the compactified picture now, we consider the equation 

f-A,+±J^G = 47ny(A), (5) 

where 5(A) is the (^-distribution centered at A, and g ab is the regular metric 
g ab = u) 2 h ab on E, the one-point compactification of E,. 

The function G is required to be positive and smooth on E\A. Near A such 
a solution will blow up, in particular it is known || that 

G = uj- 1/2 + j + 0{uo l/2 ), m = const., (6) 

where u is a RADF near A and the constant m is independent of the choice of 
RADF. 

The manifold (E, h) with h ab = uj~ 2 g ab is then asymptotically flat near A with 
zero ADM mass. The metric h ab = ip 4 h ab = G A g ab however, due to Eq. (|6|), is 
asymptotic to the Schwarzschild metric of mass m: 

^=^/2 G= l + ^l/2 + OM; (7) 

and has zero scalar curvature since (— + ^Rg)G = on E. Its mass m is 
guaranteed to be positive by the positive mass theorem |7]]. 

2.2 The General Case 

We will now incorporate extrinsic curvature into the formalism. For vanishing 
trace of the extrinsic curvature (i.e. a maximal slice) the momentum constraint 
D a K ab = is conformally invariant if K ab is rescaled to K ab = ifj 10 K ab . This 
holds in particular if the conformal rescaling corresponds to compactification. 
Here we start with a symmetric tensor P ab on the compact base manifold, which 
is transverse and trace-free (TT) with respect to the compactified base metric g ab . 
Then we obtain the TT tensor used in the AF setting by rescaling to K ab = uP ab , 
which enters the Hamiltonian constraint in its usual formulation 

f-A- h + ^R-^j = ^K ab K ab r\ i> - 1 at oo. (8) 
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The physical extrinsic curvature is then obtained by the rescaling K ab = G~ P a b- 
In order to compute the asymptotic behavior of P a t, we have to take into account 
both the rescaling, and the fact that the tensor components will be evaluated in 
different coordinate systems, which yields 

= = ™™™ G -2 = p fX (u 3 ), (9) 

3 dx l dx3 3 dx l dx3 3 13 v h y ' 

where Kij are components in an asymptotically Cartesian coordinate system x l , 
and Pfj are components in a coordinate system x 1 which is regular on S 3 near A. 
For asymptotically Schwarzschildian data (K ab = 0(lu 3 ^ 2 )) we therefore choose 
Pab = 0(lu~ 3 / 2 ). In the more general case allowing for linear momentum we have 
P ab = 0(lu- 2 ). 

In order to generalize the Hamiltonian constraint to the non-time-symmetric 
case, we consider the equation 

A- g +^R s j H = 4tt5(A) + H- ll -P ab P ah . (10) 

This can be treated as follows: denote by G a solution to the time-symmetric 
Hamiltonian constraint and define 

H = G + (f). 

This yields the following equation for the function 

L- g <p = 1(1 + G-^)- 7 G- 7 P ab P ab . (11) 

o 

Note that this equation is regular for Schwarzschildian initial data, since then 
G~ 7 P ab P ab = Oiu) 1 / 2 ). It follows from the Appendix of f| that there exists a 
unique positive solution H of Eq. ([□]), smooth on £ — A, bounded and with 
bounded first derivatives on M. 

Existence and uniqueness of solutions to the Hamiltonian constraint for given 
[Pab, 9ab) can be inferred from the time-symmetric case. It is known that a unique 
positive solution for (|5|) exists, iff the Yamabe number of g ab is positive, or, 
equivalently if the lowest eigenvalue of Lg is positive. This result was effectively 
proven by Cantor and Brill [TD|. Existence can thus be read off from the solution 



of the linear equation 

Lg0 = \6, (12) 
and is in particular independent of the choices of P ab and A. 
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3 Toroidal Conformal Symmetry and the Brill 
Ansatz 



In this section we will study metrics on S 3 in the context of toroidal symmetry. 
To start we write the standard metric on S 3 , the one-point compactification of 
R 3 , in toroidal coordinates (see e.g. Jill) as 



( /,- = + (1 - p 2 ) dtp 2 + p 2 dz\ (13) 

1 — p z 



where < p < 1, < z, <p < 2n. The vector fields (d/dip) a , (d/dz) a form a pair 
of commuting, orthogonal and hypersurface-orthogonal Killing vector fields of g 
spanning the surfaces of constant p. For p ^ 0, 1 these are flat tori, as can be 
seen immediately by inspection of the metric. Decompactifying to the flat metric 
by stereographic projection a torus p = const, is mapped to a standard (non-flat) 
torus, for <p = const.. The sets p = (respectively p = 1) are linked great circles 
on S 3 , corresponding to the z-axis (respectively the circle z — 0, r — 1) after 
stereographic projection. For more details see e.g. fill , || or 0. 

In order to achieve a convenient form of the Hamiltonian constraint for axially 
symmetric initial data, we use the well known Brill ansatz [Jl~2| |, which is usually 



formulated in the asymptotically flat picture, but straightforwardly translates to 
h = ^(c 2a ^(dp 2 + dz 2 )+p 2 d V 2 ) 



G A (e 2 ^ + (1 - p 2 ) dz 2 ) + p 2 V) = ^- 2 g- 



(14) 



Here we have chosen the conformal factor for compactification as u = Q = 
2/(1 + r 2 ), which corresponds to stereographic projection, and is compatible with 
the symmetries we want to exploit. Note that Q is actually not an RADF, but 
2Q is. For the base metric g we can now read off the ansatz 



- = part + _ 2) d _ 2 ^ + _ 2d ^ 



(15) 



The crucial step is to choose q = q(p). This implies that {d/dz) a is a Killing 
vector field (KVF) of g a b in addition to axial symmetry. At the axis p = the 
two conditions g(0) = 0, q'(0) = are required to ensure both regularity on the 
axis and the correct asymptotic behavior (see [0 and ||). A deformation of 



the standard metric on S* 3 with q = q(p) now respects both KVF's (<9/<9<p) a and 
(d/dz) a . The scalar curvature simplifies to 

f (p)=e- 2a «( 3 - + V(p)) (16) 
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where V(p) = — | [(1 — p 2 ) q" — 2pq']. The action of the conformal Laplacian on 
axially symmetric functions is thus given by 

L- g f(p, z) = {- A + 1+Vjf. (17) 

This simplification will be the key to turn PDEs into ODEs in the next section. 



4 The Constraints in the Light of Conformal 
Symmetry 

4.1 Existence of Solutions to The Hamiltonian Constraint 

Following Sec. |2.2| , to determine the existence of a solution to the Hamiltonian 
constraint is equivalent to positivity of the lowest eigenvalue Ai of L s . Assuming 
toroidal symmetry, the eigenvalue problem 

Lg6 = X6, 

is separable (since q = q(p), hence the potential depends on p only) and the 
solutions can be written as 

9\mn(p, <f, z) = &\mn(p) sin(rup + a) sm(mz + (3), 
where a and (3 are arbitrary angles. We are thus left to solve an equation for 

F 1 — p l p l J 

where Lp is the radial part of the conformal Laplacian Lg, 

2x d 2 l-3p 2 d 3 rr , . 

L P=-v-f)op — j^Tp + l + v{p) - (18) 



Since the first eigenfunction is non-degenerate |13|| and has therefore to depend 



trivially on the angular coordinates ip and z, we may set m = n = when looking 
for the lowest eigenvalue Ai : 

Lp&xaa = A^aoo- (19) 

The existence question has thus been reduced to an ODE of Sturm-Liouville 
type which can be solved by standard methods. 
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4.2 The Time- Symmetric Hamiltonian Constraint 

We will now use the setup developed in Sec. (|j) to convert the time-symmetric 
Hamiltonian constraint into a set of uncoupled ODEs. Using q(A) = 0, so that 
exp(2aq)5(\) = 5(A), the constraint equation reads 

(- A + \ + v) G(j>, z) = 4tt y/25 (A) , (20) 

where the factor of \[2 was introduced to compensate for Q not being an ADF, 
so that the solution of the unperturbed (a = 0) case becomes G = fi~ 1//2 . 
Nevertheless we will sloppily call G a Green function. 

Before we actually transform fl2~0|) to ODEs, we will formulate it as an equa- 
tion for a regular function instead of being an equation for the singular G, whose 
asymptotic behavior near A was given in (^J). To remove the distributional charac- 
ter of this equation we write the solution as a perturbation of the Green function 
of the conformally flat metric: 

G = G + (f) (21) 
and get an equation for the continuous function (p: 

-A + l + v\<p = -VG . (22) 

Depending on the choice of q the inhomogeneous term VGq still may show a 
singular behavior. If we have q = 0(p m ) then V = 0(p m ~ 2 ) and 

VG = O( 1 



'z* + p 1 



which is singular for m < 3. Later we will cosine transform with respect to z and 
thus have to consider finiteness of the integral 

VG dz 

which will be finite for m > 2 (Vo := V"(0) = 0). The boundary condition 
q = 0(p 2 ) translates to m > 2, so m = 2 is both the minimal value allowed 
by the boundary condition and the critical value for the finiteness of the Fourier 
transform of the inhomogeneous term. Note that for the critical value the diver- 
gence is logarithmic. To further regularize the equation for m = 2 we make the 
ansatz 

g = ir 1/2 + Cl n 1/2 + c 2 n 3/2 + 0. (23) 
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Now the constants Cj, c 2 are chosen so that the contributions of the fl 1//2 and 
fi 1 / 2 terms cancel at A. The necessary choices are 

ci = V , c 2 = \v (2 + V ). 
o 

The equation for then becomes 

, 2 , d 2 l-3p 2 d 2 1 d 2 3 *\ j , 

-«^ + ^ap + 1 "^J* = < 24 > 

where the function F is given by 

f = (v — v ) n- 1 / 2 + v {v - v ) n 1 ' 2 + ^ (2 + v ) (6 + v) ft 3/2 , (25) 

o 

which reduces to 

F = Vtt~ 1/2 (26) 

for q = 0(p 2+m ). 

The simple form of fl24j) , namely that the conformal Laplacian does not explic- 
itly depend on z, suggests the method of expanding <f>(p, z) into an even Fourier 
(cosine) series {(f) is invariant under reflection at the equator) 

00 

(j)(p, z) = ^ ^nip) cos nz. (27) 

n=0 

Applying an analogous transformation to F gives a sequence of uncoupled linear 
ordinary differential equations for the Fourier coefficient functions <p n : 

P " « W + ^1 " T^V " i " K ) *•<*> = W (28) 



Equation ( gSH has singular points at p — 0,1, which can be associated with the 
two axes of symmetry and the corresponding coordinate singularities. Boundary 
conditions are necessary to guarantee that the solutions are regular. 

The axis at p = is also a physical axis of symmetry, and the boundary 
conditions there directly correspond to the ones in the asymptotically flat context. 
In addition to requiring regularity of q and the conformal factor ip (q,if) G C 2 (IR 2 ), 
this implies q p and ip p both have to vanish on the axis. From ( p3|) and the 
definition of ip as ip = VL^^G (where fl is the compactifying conformal factor) we 
get 

1>(jj(x, y, z), z(x, y,z)) = l+V n+^-(2 + V ) Q 2 + fi 1/2 

o 

so that instead of working with ip we impose equivalent conditions on <fi(p(x, y, z), 
z(x,y, z)), which is simply the sum of the Fourier terms T n = <f) n (p) cosnz. For 
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p < 1 the coordinate transformation from cylindrical coordinates to toroidal 
coordinates is regular, so we simply require 4> n (f>) to be a regular function in this 
interval with 0' = at p = 0. 

By contrast the line p = 1 is only an axis of symmetry for the compact 
manifold and the coordinate singularity there does not bear any physical meaning, 
and indeed the situation at p = 1 is more delicate. The coordinate transformation 
from (p, z) to (p, z) is singular there. The boundary conditions on the functions 
4> n {p) at p = 1 have to ensure that the nth Fourier term T n = 4> n {p) cos nz is a 
regular function of the coordinates x, y, z. One can show that (see M) a necessary 
and sufficient condition for regularity of T n is that 

/ — \ Mp) 



is a regular function. 

We also introduce, analogously to <p n , the functions 



-n 



We get 



i, 2\ d 2 1 - (3 + 2n)p 2 d 2 n 3 _ \ . , on . 

The boundary value problem of finding a positive Green function has thus 
been translated into the problem to find C 2 solutions to ( |29|) for every n. An 
approximate solution is constructed by solving for n < N, resulting in 



G N = IT 1/2 + y ^ 1/2 + y (2 + Vb) fi 3/2 + ^ v 7 ! - P 2 V cos n^. (30) 

n=0 

An efficient method for solving one-dimensional boundary value problems is the 
shooting and matching technique, which has been used in the numerical study 
discussed in 0. 

So far we have only explicitly dealt with the axially symmetric case, by assum- 
ing that the ^-distribution source was placed on the axis p = 0, in particular at 
z = 0. This renders G and <p axially symmetric and in the Fourier series we only 
have to consider the cos nz terms. The generalization of the axially symmetric to 
the general 3D case is straightforward. The only places where axial symmetry has 
been used so far is in the placing of the source, and the corresponding functional 
form of the Green function of the conformally flat conformal Laplacian G (i, A), 
where A is placed on the axis p = 0, which is 



N 



Gq = (1 — a/ 1 — p 2 cos z) 
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For a general position of the source at (p', z', p') we instead have that 

Go = (1 - V1-PV1-P' 2 cos (z - z') - pp' cos (<p - cp')). 

This can again be Fourier expanded in a double Fourier series in z and ip. 

Note that the structure of KVF's of the physical metric depends on the choice 
of "point at infinity" . The condition has been formulated as a theorem by Beig 
[fPH : A KVF rf of a metric g a b is also a KVF of the decompactified metric G 4 g a b, 
where G satisfies the time symmetric Hamiltonian constraint, L g G = Att5A, with 
A the point at infinity if and only if A is a fixed point, i.e. r] a \x = 0. 

Since the axes will turn out to be disjoint in the present case, this means that 
only one, or even none, of the KVF's may survive decompactification. 



5 Initial Data Containing Black Holes 

The formalism of solving the constraints on a conformally compactified back- 
ground discussed in Sec. |2| and the simplifications due to assuming conformal 
symmetry described in Sec. [| can be immediately generalized to the case when 
there are several asymptotic ends. We simply choose a finite number of points 
A, (z = 1, . . . , N), solve numerically for the corresponding Green function Gi 

Lgd = 47TCi5{Ai), 

and define 

hob = I y \ c i G i ) g abj (31) 



where the "source strengths" Cj are arbitrary positive numbers. When the am- 
plitude parameter a vanishes, these are simply the many-black-hole solutions 
discussed by Brill and Lindquist [jT5|. The existence of the Gi is guaranteed by 



a positive Yamabe number of g. Thus, starting from regular initial data (with 
positive Yamabe number), we can always add as many black holes as we like, 
with arbitrary weights Cj and positions Aj. While the Green function yields the 
solution to the time-symmetric constraint, the generalization to incorporating ex- 
trinsic curvature is again straightforward along the lines of Sec. [2.2| , in particular 
one ends up with Eq. (|TT|) . A similar approach to treat black- hole initial data, 
which does not use conformal compactification, has recently been discussed by 



Brandt and Brugmann |16j in the conformally flat case. 

If all the Aj lie on the two axes, each Gi in G = Yli=i c i^i corresponds to an 
axially symmetric (Brill-type) solution. If at least one of the sources is placed 
off-axis, the resulting physical geometry will not have any Killing vector fields at 
all. 

Starting with a globally regular Brill-wave solution and adding sources, it 
is interesting to observe that Brill's simple positive-mass proof ]I2| immediately 
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gives positivity of mass also in this non-axially symmetric, topologically nontrivial 
situation: For a single source placed at Ai we have that 

G 1 = G + ^- + O(G l ), 



where positivity of m\ is guaranteed by Brill's argument ||12|| . Adding sources we 
get a different mass fhi at Ai, 



N 



m l = 2{G-G )\A 1 =m 1 + 2j2Gi\Ai 



n=2 

where the second term is positive by positivity of each individual G{. 

Although d/dz is not a KVF of the physical metric hab = G 4 g a b, i.e. arbitrary 
rotations in z are not symmetry operations, h a b will possess a discrete symmetry 
of invariance under coordinate transformations z — > z + 5z if the Aj are placed 
equally spaced at 

Zj = , Sz = (32) 

1 N+V N+l K 1 

and are all equally weighted, that is Wi — 1 with the choice wq = 1 or more 
generally for any arrangement of sources that is periodic in z. 

This can be considered a generalization of inversion symmetry for two asymp- 
totically flat sheets which corresponds to placing sources of equal strength at 



antipodal positions. Like any symmetry of the initial data [17] this symmetry 
is respected by time evolution. A nice consequence is that the symmetry carries 
over to minimal surfaces, that is if one finds one minimal surface, one can get 
another N —1 minimal surfaces by shifting the coordinate z. In particular every 
asymptotically flat region will show an apparent horizon. This means that for 
every choice of the function q the string of black holes will be surrounded by a 
common apparent horizon. 

A particular useful formulation of toroidal symmetry is that, in the case of a 
general distribution of Aj's, not all of the G^s have to be computed separately: 
for any A 1; A 2 lying on the same orbit of the isometry group, G2 is given in terms 
of G\ by a product of suitable rotations in (p and z. Take e.g. N = 2 with Ai at 
(p = 0, z = 0) and A 2 at (p = 0, z — a). 

In the special case a = tt, the physical metric h a b has the mirror symmetry 
z 1 — ► 2tt — z. When c\ = c 2 , it has the inversion symmetry {z} U {z + 7r} 1— > 
{a — z} U {a — z + 7r}, leaving the "throat" {z = a/2} U {a/2 + tt} invariant. 
Thus the "throat" is totally geodesic with respect to h a t, in particular a minimal 
embedded 2-sphere (a "horizon"). When h a b is conformally flat, i.e. a = 0, 
the above discrete symmetries are present irrespective of ci, c 2 and a. Because, 
then, by combining a homothety with a proper conformal motion coming from 
the Killing vector on (i? 3 , 5 a6 ), we can always arrange for c[ = c' 2 , a' = n. The 
corresponding physical data is of course nothing but a time-symmetric slice of 
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the Kruskal spacetime with mass m(ci, 02, oi) = 2v^2ciC2(l — cosct) -1 / 2 . When 
a ^ our solutions for c± 7^ C2 are not inversion symmetric. Placing sources of 
equal strength at p = 0, a>i = ^r, i = 0, 1, . . . , N — 1, we find in an analogous 
manner that on R 3 \ IJILi ^% eacn asymptotic end is surrounded by a throat. 
Thus, viewing Ai as "infinity" and the other Aj's as "black holes" or "particles" 
we can in particular say that these objects, viewed from infinity, are so close 
together that, in addition to a horizon for each of them individually, there is one 
surrounding them all. 
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